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SECTION  1 

INTRODUCTION 


The  term  adaptive  optics  generally  refers  to  the  use  of  deform¬ 
able  mirrors  to  control  phase  aberrations  that  are  involved  in  the 
generation  of  electromagnetic  radiation  and  its  propagation  through 
a  medium  such  as  the  earth's  atmosphere.  Such  control  is  important 
both  for  passive  celestial  observation  through  the  atmosphere  and  for 
laser-beam  propagation,  especially  high  power.  The  latter  application 
is  the  main  motivation  for  this  report. 

There  are  several  physical  phenomena  that  give  rise  to  the  need 
for  adaptive  optics,  including  beam  quality,  jitter,  turbulence,  and 
thermal  blooming.  Of  these,  all  but  the  last  are  said  to  be  linear: 
correcting  the  phase  aberrations  they  introduce  does  not  in  any  way 
feed  back  to  change  the  phenomena  causing  the  aberrations.  Thermal 
blooming  is  different:  it  consists  of  a  lensing  effect  in  the  atmosphere 
caused  by  the  heating  induced  by  the  propagating  beam,  which  in  turn 
is  modified  by  any  phases  introduced  at  the  source. 

Much  work  on  theoretical  formulation,  computer  simulation,  and 

experimental  demonstration  has  been  undertaken  in  order  to  arrive  at  a 

(1  2  3 )  * 

workable  adaptive  optics  capability.  '  '  '  So  far,  results  are 
rather  less  than  what  might  have  been  hoped  for.  Apparently,  diffrac¬ 
tion  limited  performance  cannot  be  guaranteed  in  the  presence  of 
thermal  blooming. 

Given  the  suspected  disappointing  behavior,  there  has  been  a 
pressing  need  for  a  physical  understanding  of  the  fundamental  limita¬ 
tions  involved  in  the  adaptive  optics  concept.  Furthermore,  the 
engineering  procedures  for  implementing  adaptive  optics  need  to  be 
thoroughly  reviewed  in  the  light  of  this  new  understanding.  This 


* 

Superscript  numerals  refer  to  similarly  numbered  items  in  the  List 
of  References. 
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report  attempts  to  provide  the  necessary  analysis  at  a  level  simple 
enough  to  facilitate  physical  understanding,  but  accurate  enough  to 
assist  in  contriving  remedies  to  the  extent  possible. 

The  analysis  is  based  upon  an  algebraic  model  for  continuous- 

wave  (CW)  thermal  blooming  reported  earlier  by  Breaux,  Evers,  Sepucha, 

(4 ) 

and  Whitney.  This  model  was  developed  for  the  purpose  of  providing 
a  standard  hand  calculation  tool  for  systems  analysis  exercises,  where 
the  numerous  parameters  are  to  be  explored,  making  more  detailed  computer 
simulation  of  atmospheric  propagation  impractical.  Its  key  feature  in 
regard  to  thermal  blooming  is  characterization  of  blooming  spread  in 
terms  of  an  analytic  integral,  which  represents  the  accumulation  along 
the  range  (starting  from  the  center  of  the  aperture)  of  phase  perturbation 
due  to  heating  in  a  beam  with  absorption,  scattering,  convective  clear¬ 
ing,  and  focusing. 

The  original  algebraic  model  required  several  extensions  for 
application  to  the  adaptive  optics  problem.  First,  the  original  model 
assumed  round,  nearly  Gaussian  beams  on  target,  and  characterized  all 
beam  spreading  in  terms  of  a  single  1/e  radial  dispersion.  In  reality, 
blooming  produces  a  beam  which  is  not  round,  and  in  fact  has  only  one 
axis  of  symmetry  (a  sort  of  "light-bulb"  shape).  This  means  that  a 
number  of  Zernike  modes  are  involved.  The  present  model  has  been 
extended  to  include  appropriate  modes  up  through  spherical  aberration. 
Second,  the  original  algebraic  model  did  not  include  deterministic 
phase  aberrations  introduced  at  the  aperture  to  attempt  compensation  of 
the  blooming.  The  present  model  includes  these,  plus  the  resulting 
change  in  the  blooming  spread  caused  by  its  nonlinear  dependence  on 
what  is  applied  at  the  aperture.  Third,  the  original  model  was  purely 
static.  The  present  model  allows  for  temporal  variation  of  the  control 
variables. 

Section  2  extends  the  static  model,  and  Section  3  notes  the  impli¬ 
cations  evident  at  that  point.  The  major  implication  available  is  a 
clear  statement  of  the  fundamental  limits  on  adaptive  correctior.  of 
thermal  blooming.  The  numbers  are  not  impressive:  a  factor  of  2.6  for 
uniform  round  beams,  improvable  only  to  a  factor  of  8  by  manipulating 
the  beam  shape.  These  results  mean  that  any  investment  in  the  develop¬ 
ment  of  adaptive  optics  for  CW  thermal  blooming  must  be  motivated  not 
by  the  hope  of  correcting  that  problem,  but  rather  by  the  hope  of  reduc¬ 
ing  its  interference  with  correction  of  the  simpler  linear  phenomena. 


Even  though  the  payoff  for  optimal  thermal  blooming  is  modest, 
the  penalty  for  doing  it  conventionally  is  potentially  severe.  Section  4 
shows  that,  if  applied  to  thermal  blooming,  any  of  the  classic  adaptive 
optics  concepts  commonly  considered  in  the  past  could  be  counter  produc¬ 
tive  and  end  up  causing  a  net  degradation.  The  section  also  shows  how 
the  destructive  behavior  could  be  remedied. 

The  thermal-blooming  compensation  techniques  implied  by  the  static 
analysis  of  the  phenomenon  are,  unfortunately,  inappropriate  for  the 
linear  phenomena.  The  requirements  of  linear  and  nonlinear  phenomena 
tend  to  conflict.  The  possibility  for  separating  the  phenomena  so  that 
each  may  be  dealt  with  appropriately  resides  in  dynamics.  Section  5 
introduces  dynamic  models  and  explores  their  implications. 

Section  6  concludes  with  recommendations  for  verifying  the  findings 
of  this  report. 
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SECTION  2 


STATIC  THERMAL-BLOOMING  ANALYSIS 


The  static  analysis  constitutes  a  propagation  model,  showing  the 
effect  of  aberrations  upon  peak  intensity  and  discussing  their  physical 
generation . 

2 . 1  Mathematical  Model  for  Aberrated  Beams 

To  reach  the  heart  of  the  thermal-blooming  problem  as  it  relates 
to  adaptive  optics,  we  require  a  model  with  reasonable  fidelity,  facility 
to  vary  applied  phase,  and  analytic  tractability .  Models  available  to 
choose  from  include  complex  wave-optics  computer  codes  that  produce  a  beam 
profile  in  the  target  plane,  more  simplified  codes,  and  various  analytic 
approximations,  usually  dealing  only  with  peak  intensity.  The  require¬ 
ments  are  best  met  by  a  simple  algebraic  model  along  the  lines  of  Ref¬ 
erences  4  and  5.  This  model  represents  the  results  of  efforts  to  fit 
scaling  laws  to  a  large  body  of  data  generated  by  detailed  atmospheric 
propagation  simulations  done  at  the  Naval  Research  Laboratory  (NRL) 
using  a  finite-difference  wave-optics  code. 

Any  such  wave-optics  code  provides  a  detailed  intensity  profile 
on  target,  from  which  several  different  summary  characteristics  can  be 
extracted.  These  include,  for  instance,  peak  intensity,  line-of-sight 
beam  dispersion,  and  beam  area.  The  last  of  these  further  admits 
several  definitions,  including  area  to  some  percentage  of  total  power, 
and  area  defined  by  a  ratio  of  squared  integral  of  intensity  to  integral 
of  intensity  squared  (suggested  independently  by  Lincoln  Laboratory  and 
Draper  Laboratory  researchers).  If  beams  on  target  were  Gaussian  in 
shape,  all  such  characteristics  would  convey  equivalent  information. 

The  simplified  algebraic  model  assumes  that  this  is  nearly  the  case, 
and  speaks  nominally  of  peak  intensity  on  target. 


as 


The  original  algebraic  model  approximates  peak  intensity  on  target 


Ip  =  PT/ttZ2  o2 


2 

where  P  is  power,  T  is  transmission,  Zfc  is  target  range,  and  o  is  the 
1/e  radial  beam  spread  due  to  the  combination  of  all  dispersive  effects, 
including  diffraction  and  beam  quality  (D) ,  high-frequency  turbulence 
(T) ,  high-frequency  jitter  (J) ,  low-frequency  wander  (W) ,  and  blooming 
(B). 


a 


2 


+ 


a 


2 

B 


The  above  model  applies  for  beams  that  are  nearly  round  on  target. 
More  generally,  we  can  allow  an  arbitrarily  shaped  beam  with  principal 
axes  x  and  y,  so  that 


In  addition,  the  model  above  contains  no  deterministic  ray-deviating 
effects  other  than  thermal  blooming.  More  generally,  we  can  allow 
arbitrary  aberrations  applied  at  the  aperture  by  the  replacement 


a 


2 

B 


a 


2 

R 


where  R  refers  to  the  net  result  of  all  deterministic  ray-bending  effects. 

Since  aberrations  at  the  aperture  are,  literally  speaking,  "linear"  effects, 

2 

we  introduce  the  notation  o^d)  to  refer  just  to  the  random  linear  effects: 


°L(1> 


2  2  2 
°D  +  °T  +  °J  + 


o 


2 

W 


Taking  account  of  the  dissimilar  principal  axes, 


2  2  2 

o  =  0,  (1)  +  0„ 

X  Lx  Rx 


2  2  2 

a  =  a,  (1)  +  at 
y  Ly  Ry 


5 


The  dispersions  due  to  ray-deviating  effects  depend  on  whatever 
Zernike  mode  aberrations  are  present.  A  general  formulation  can  be 
built  up,  after  first  considering  some  special  cases,  corresponding 
to  each  type  of  aberration  that  blooming  might  present  (see  Table  1) . 
For  each  case,  we  shall  draw  a  geometric  ray  model  (see  Figure  1),  and 
formulate  an  algebraic  spot  model.  Then,  we  shall  formulate  a  general 
ray  model  and  spot  model  covering  all  aberrations  and  combinations 
thereof. 


Table  1.  Zernike  optical  modes. 


i 

P^r^) 

Name 

0 

1 

1 

r  cos  9 

x-tilt 

2 

r  sin  9 

y-tilt 

3 

2r2  -  1 

defocus 

4 

r2  cos  29 

astigmatism 

5 

r2  sin  29 

astigmatism 

6 

(3r3  -  2r)  cos  9 

primary  coma 

7 

(3r3  -  2r)  sin  9 

primary  coma 

8 

r3  cos  39 

trefoil  coma 

9 

r3  sin  39 

trefoil  coma 

10 

6r4  -  6r2  +  1 

spherical 

Notes:  1.  Variable  r  is  normalized  to  unity  at  the 

aperture  radius. 

2.  Measuring  9  from  the  direction  of  wind,  the 
modes  of  interest  are  3,  4,  6,  10. 

Let  us  begin  with  Focusf6^  Allowing  a  constant  m"  to  account 
for  different  beam  shapes/5^  we  have 

°Fx  '  4y  ' 

where  NR  is  the  ratio  of  target  range  to  focal  range 

NR  =  Zt/Zf 
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FOCUS: 

RAYS  CONVERGE  TO  A  POINT  OFF  TARGET 


ENLARGEMENT 


ASTIGMATISM: 

TWO  FOCAL  LENGTHS  EXIST.  ONE  FOR  x,  ONE  FOR  y 


^ELONGATION  PERPENDICULAR 
TO  WIND 


COMA: 

FOCAL  LENGTH  IS  A  FUNCTION  OF  x 


BENDING  INTO  THE  WIND  DUE 
TO  RAY  CROSSING 


SPHERICAL  ABERRATION: 

FOCAL  LENGTH  IS  A  FUNCTION  OF  RADIUS 


ENLARGEMENT 


Z.  z. 


Figure  1.  Ray  pictures  of  Zernike 
mode  aberrations. 
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In  the  case  of  astigmatism,  there  are  two  focal  ranges  instead 
of  one,  and  we  have 

4  - 

4y  *  <0-S»-  VZt>2(1-%'2 

where 

"ft.  ■  V2f* 

NRy  =  Zty/Zfy 

When  the  astigmatism  is  caused  by  blooming  of  a  beam  nominally 
focused  on  target,  it  will  be  accompanied  by  an  overall  defocus  as  well, 
so  that  in  effect  there  is  an  elongation  of  the  beam  perpendicular  to  the 
wind,  as  noted  in  Figure  1. 

Coma  is  characterized  by  having  one  axis  perturbed  in  such  a  way 
that  the  rays  never  come  to  a  point.  There  is  a  perturbation,  but  much 
smaller  for  the  other  axis.  Therefore,  coma  is  well  represented  by 

4  *  ‘°-5m-  W 2(1 -"a1 2 

4  -  ..sfl-K,.)2*!!-!!,,)2] 

where 

“a-  '  V2- 
V  '  V2* 

When  the  coma  is  caused  by  blooming  of  a  beam  nominally  focused 
on  target,  it  too  will  be  accompanied  by  an  overall  defocus.  This 
means  that,  in  the  target  plane,  most  ray  crossing  is  displaced  slightly 
into  the  wind,  causing  the  well  known  "bending  into  the  wind"  associated 
with  blooming  (see  Figure  1) . 


Spherical  aberration  is  axially  symmetric,  but  the  rays  never  come 
to  a  point  so  both  factors  in  the  dispersion  are  never  zero: 


4  ■  4 


(0.5m"Rni/Zt)2(0.5)  [(1-NR_)2  +  (1 


-W2] 


where 


nr-  =  Vz- 


R+ 


-  vz+ 


Consequently,  spherical  aberration  always  enlarges  the  spot  delivered. 

The  dispersion  formulas  for  the  various  aberrations  previously 
offered  do  follow  a  general  pattern.  Regardless  of  the  aberration,  we 
can  model 


4  -  0-5<4+4l 


4y  ■  0-5"4  +  »R4* 


where 


Ri 


(0-5”-  W  2 11 -Si’2 


with 


-  Vzfi 


The  particular  cases  are  recovered  by  the  substitutions  noted  in  Table  2. 

Normally,  the  parameterization  of  an  aberrated  beam  is  in  terms 
of  Zernike  coefficients.  The  aberration  phase  y  experienced  by  a  ray 
originating  at  the  point  r,0  in  the  aperture  is 


’l'<r,e)  =  k  £  Cz.  P.  (r , 6 ) 

i 
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Table  2.  Special 

cases  of 

general 

ray  model. 

FOCUS 

Zf  - 

Z1 

=  Z2  = 

Z3  = 

Z4 

ASTIGMATISM 

(Z1  = 

Z2  > 

¥  (z3 

II 

N 

COMA 

(Zf  = 

Z1 

=  Z2> 

¥  (z3 

¥  z4) 

SPHERICAL  ABERRATION 

(Z!  =  z3>  *  (Z2 

= 

where  k  is  wave  number  2ir/X  (with  X  being  wavelength)  ,  and  C„ .  is  the 

t  h  ®  - 

coefficient  of  the  i  Zernike  mode.  The  Zernike  coefficients  of 

interest  in  this  model  are  in  one-to-one  correspondence  with  linear 

combinations  of  the  dispersion  parameters.  The  relationships  are  as 

follows: 

(1)  Because  an  aberration  in  focus  means  all  rays  are  similarly 
converged  off  target,  the  corresponding  Zernike  coefficient 
Cz3  is  proportional  to  the  average 

(0Rl  +  0R2+aR3+aR4>/4 


(2)  Because  astigmatism  means  a  difference  between  x  and  y  axes, 
the  corresponding  Zernike  coefficient  Cz4  is  proportional  to 


(°R1  + 


R2 


-  a 


R3 


-or4)/4 


(3)  Coma  implies  a  dissimilarity  between  maximum  and  minimum 
focal  ranges,  but  only  on  one  axis,  not  the  other.  The 
corresponding  Zernike  coefficient  Cz6  is  proportional  to 


(0Rl-0R2-°R3  +  aR4)/4 
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(4) 


W*" — l'^»i  ■  . . 

»-w- 


Spherical  aberration  is  an  axially  symmetric  dissimilarity 
between  maximum  and  minimum  focal  ranges.  The  corresponding 
Zernike  coefficient  CZ1Q  is  proportional  to 

(aRl-°R2  +  aR3-aR4)/4 

For  each  of  these  modes,  the  actual  proportionality  constants 
are  functions  of  beam  shape.  The  simplest  example  is  the  uniform  beam. 
For  this  case.  Reference  7  suggests  that  the  proportionality  is 


ai  =  °Dk(  /Pir  dr)1/2 

Notice  the  behavior  of  the  model  in  the  limit  of  axially  sym¬ 
metric  random  effects  and  of  very  small  aberrations.  For  that  case, 
the  beam  dispersion  is  equivalently 

«2  -  O2(1)*0.5(»21*^2*«2,h.<,24I 

=  oj(l)  +0.5(Op+o^+a^+Og) 

because  the  various  aR  cross  terms  in  the  squared  Zernike  a's  cancel 
each  other.  That  is,  whether  viewed  in  terms  of  dispersions  or  Zernike 
coefficients,  the  model  parameters  do  not  cross  couple. 


2 . 2  Interaction  of  Ray-Bending  Effects 


The  ray  deviation  effects  represented  by  the  dispersions  o  .  or 
the  Zernike  coefficients  Czi  are  the  net  result  of  aberrations  imposed 
at  the  aperture  (F)  and  aberrations  induced  in  the  atmosphere  by 
thermal  blooming  (B) .  The  parameters  describing  the  blooming  are 
complicated  nonlinear  functions  of  the  parameters  of  the  aperture 
aberrations  (and  really  even  each  other) .  But  the  heart  of  the  adaptive 
compensation  concept  is  that  despite  nonlinearities,  aberrations  applied 
at  the  aperture  must  to  some  extent  cancel  those  in  the  atmosphere.  In 
terms  of  Zernike  modes,  one  might  suppose 


'Zi 


CFi  +  CBi 
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The  concluding  remarks  of  Section  2.1  show  that  in  the  limit  of  small 
total  aberrations,  this  would  be  equivalent  to 


Ri 


°Fi  + 


Bi 


In  an  operational  adaptive  optics  system,  the  blooming  parameters 
would  in  some  sense  be  measured  (say,  in  a  Zernike  mode  expansion  of  a 
measured  phase  profile) .  But  for  the  purposes  of  this  analysis,  which 
intends  to  anticipate  the  behavior  of  such  systems,  a  mathematical  model 
is  needed.  It  should  be  understood  that  what  such  a  model  lacks  in 
exactness  can  be  remedied  by  measurement  in  a  real  system. 

We  begin  with  the  scaling  law  formulations  found  in  References  4 
and  5,  which  model  the  total  beam  dispersion  in  the  absence  of  aperture 
corrections  as 


2  2  2 
a2  =  o2(l)+a2 

Compare  this  to  the  present  model  evaluated  in  this  limit 


2ax  ay 


2  [  (0 . 5a2  (1)  +a2x)  (0.5o2(l)  ^  ^)}1/2 


Numerical  similarity  of  the  two  models  requires,  approximately 


2  2 
°Bx  ~  °By 


o2/2 


The  above  approximate  equality  between  crBx  and  oBy  is  further 
confirmed  by  a  ray  analysis  following  Figure  1.  Figure  2  shows  stylis¬ 
tically  the  rays  at  x  *  ±Rm'  y  =  ±Rtn  as  deviated  by  the  blooming.  The 
ray  at  x  =  -Rm  is  undeviated  because  of  wind  clearing.  The  rays  at 
y  =  ±Rm  are  somewhat  deviated.  The  ray  at  x  =  +RJn  is  twice  as  deviated 
due  to  greater  heating.  The  geometry  suggests  the  approximate  equality. 
In  addition,  it  suggests 


°B1 


0 
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y 
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Figure  2.  Rays  deviated  by  bloomina. 


which  implies 


°B2  ”  °B 

and,  at  least  for  a  uniform  round  beam  shape,  it  suggests 

°B3  ~  °B4 


which  implies 

°B3  =  °B4  =  °B/>/^ 

With  the  aforementioned  geometric  arguments,  the  problem  of 

evaluating  the  blooming  parameters  is  reduced  to  that  of  evaluating  the 
2  2 

gross  o  .  The  a  can  be  represented  by  a  variety  of  scaling  law  models. 

B  B  (5) 

For  example,  there  is  the  model  of  Breaux  and  the  more  simplified 

(4) 

model  later  proposed  by  Breaux,  Evers,  Sepucha,  and  Whitney;  both 

based  on  the  phase  perturbation  due  to  heating,  measured  from  the 

aperture  center  to  target.  Both  models  are  similar  in  one  important 

2 

respect:  the  heating  phase  enters  oB  to  a  power  greater  than  unity. 

This  is  the  fundamental  feature  required  to  reproduce  the  well  known 
physical  phenomenon  of  critical  power;  that  is,  optimum  power  which 
produces  maximum  peak  intensity  on  target.  The  form  of  the  model  in 
Reference  4  is 


r'  (n 


A  V8 


Where  and  a  are  parameters  depending  only  on  beam  shape,  and  4*^ 
is  the  heating  phase  defined  in  Reference  4  as 


j. 

*h  "  /  * h 

»  rv 


(z )  dz 


where  z  =  Z/Zfc  and 


V2> 


+Hp  Np  exp(-eZ)  f dt  I  (x0  (2  .  ti  ,yQ(z,t)j 

so(RS(2,/'Vr(2,V/0w)^/l'Rm) 


with 


2.333  x  10"9  PaZ?/R3  U 
t  m  w 


Np  =  2TTRm/XZf 


nq 


(v^)1/7(VM) 


and  in  turn 


P  =  power 

a  =  absorbtion  coefficient 
Uw  =  wind  velocity 
M  =  beam  quality 

The  e  is  total  extinction,  including  scattering  as  well  as  absorption. 

The  xQ  and  yQ  are  coordinates  in  the  aperture  plane,  which  depend  on 

time  t  as  explained  in  Reference  4.  The  Rg(z)  represents  spot  radius 

at  z,  approximated  by  root  sum  squaring  ray  convergence  and  diffraction 

effects.  For  the  aberration-free  case  treated  in  Reference  4,  the  ray 

2 

convergence  contribution  has  the  very  simple  form  (1  -  z)  ,  and 


where 


(m 


2  2 
°D+°T  + 

°DOm" 


oj)1/2 


with  m  and  m"  being  constants  dependent  on  beam  shape,  and 

°D0  =  0.5(XARm)2 


The  first  step  in  generalizing  this  model  to  include  applied 

2  2 

aberrations  is  to  generalize  the  factor  aL  -  aw  everywhere  it  occurs. 
Allowing  for  everything,  we  should  have 


2  „2 
°L~  °W 


2  [ (CT2  -  02  )  (a2  -  a2  )1 

L  Lx  Wx  Ly  Wy  J 


1/2 


where 


Lx 


’U  •  °l*a)  *0.51^3*4,1 


For  scaling  laws  in  both  References  4  and  5,  the  exponent  of 

is  not  greater  than  2.  Consequently,  the  factor  in  question  enters 
2 

aB  to  a  power  greater  than  zero.  The  impact  of  the  generalization  is 
illustrated  in  Figure  3  for  the  case  where  the  applied  aberration  is 
defocus.  The  figure  shows  the  y-axis  rays  defined  by  ± (a^  +  aB3)  or 
± (aF4  +  ) .  These  are  symmetric  because  y  is  perpendicular  to  the 

wind.  Because  oB  grows  with  aF3  and  oF4»  the  rays  have  a  decreased  cone 
angle  beyond  the  focus.  This  means  the  blooming  is  acting  as  a  diverg¬ 
ing  lens,  which  is  physically  correct. 

The  next  step  in  generalizing  the  model  is  to  generalize  the 
phase  integral  The  analytic  integration  procedure  given  in 

Reference  5  applies  for  beams  focused  on  target,  and  readily  generalizes 
to  beams  focused  off  target  (ND  /  1) .  The  result  is 

K 


•  p  1 

(N_  +  N_  +  A)  (Nc  +  1) 

4-h  =  [ND  Np  Nr/N^  A]  exp(-cZt/2)  In 

-BNr  +  Ns  +  Nr  +  AC 
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OB 


Figure  3.  Gross  behavior  of  aQ. 

where 

r  2  2 1 

A  =  [(Ns+Nr)  +  (2NRNQ/NF)2j 

B  =  Ns+Nr+  (2Nq/Nf)2Nr 

1/2 

C  =  [(1-Nr)2+  (2NrNq/Nf)2] 

In  principle,  it  is  possible  to  generalize  the  expression  for 

even  further  by  using  a  more  correct  spot  model  R„(z).  What  is  required 
2 

is  that  (1-z)  be  replaced  by 

a(z)2  =  2a  (z)  a(z) 

x  y 

a2(z)  =  0 . 25 1 (1  -  Zl)2  +  (1  -  z2)2) 

a2(z)  =  0.25(  (1  -  z,)2  +  (1-z.)2) 


Where 


Such  a  substitution  can  always  be  handled  by  numerical  integration. 
However,  it  does  not  appear  likely  that  analytic  integration  would 
remain  possible  without  further  approximations. 

In  order  to  retain  an  essentially  analytical  model  we  shall  use 
the  functional  form  of  ^(Np),  and  with  linear  NR  replaced  by  the 
average  NR^f  and  quadratic  (1-NR)^  replaced  by  the  area  excess 
caused  by  the  N  .  : 

R1 


SECTION  3 
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To  solve  this  equation  formally  for  the  optimal  oF^,  let  us  expand 

to  first  order  in  o_..  Note  that  the  dependence  of  on.  on  o„.  is  such 

r  X  bl  r X 

that  is  proportional  to  : 


3  o 
To 


Bi 

Fi 


a_ .  3  a 

Fi 


Bi  3  a 


2 

Bi 

2 

Fi 


aFi  °Bi 


2  (a£  (1) 


°Fi  + 


Fj 


where  j  =  i  +  1  for  i  odd,  i  -  1  for  i  even.  The  consequence  of  this 
dependence  is  that  the  optimal  has  the  form 


-a 


Bi 


i  .3°Bi  3\  /3° Fi 

3Vh  3NRi  /  3NRi _ 

,  .  °Bi3°Bi  .  3  °Bi  34h  /3°  Fi 

aF i  3 CTF i  8Th  9NRi/3NRi 


This  expression  for  the  optimal  0°^  makes  it  clear  that  the  compensation 
usually  cannot  be  perfect.  In  addition,  it  is  interesting  to  note, 
that  with  <*0pj/3NRi  being  a  negative  number  and  everything  else  being 
positive,  it  is  entirely  possible  for  the  other  parameters  to  be 
severe  enough  that  oFi  turns  out  the  same  sign  as  oBi,  rather  than 
opposite.  This  is  a  most  startling  observation. 

Now  consider  optimization  of  °x  as  a  whole,  rather  than  the 
oRi  individually.  Because  the  formula  for  ax  is  entirely  symmetric  in 
the  various  a_.  ,  it  seems  most  likely  that  optimal  overall  correction 
ends  up  with  the  various  o„.  all  equal,  though  not  necessarily  indi- 

K 1  - 

vidually  minimum  with  respect  to  their  corresponding  In  terms  of 

Zernike  modes,  such  a  solution  would  have  all  the  high-order  modes 
corrected  exactly,  and  focus  corrected  optimally.  The  beam  would  be 
round,  and  the  rays  would  be  converged  on  or  near  target. 

Even  if  not  strictly  optimal  in  terms  of  peak  intensity,  such  a 
solution  would  be  optimum  in  more  general  terms,  such  as  its  inter¬ 
action  with  a  hot-spot  tracking  scheme.  Consider,  for  instance. 
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Figure  1,  which  shows  how  the  combination  of  coma  with  defocus  causes 
a  displacement  of  the  beam  bright  spot.  The  reverse  will  also  be  true: 
the  image  of  the  hot  spot  will  also  be  displaced  in  the  same  direction. 
Thus,  coma  causes  the  hot  spot  to  be  missensed  and  the  beam  to  be  mis¬ 
placed  with  a  net  double  error.  This  type  of  error  is  not  a  wave- 
front  tilt,  but  does  affect  the  system  like  a  tilt,  and  can  be  referred 
to  as  tilt  aliasing.  Tilt  aliasing  is  cured  in  the  optimization  des¬ 
cribed  previously. 

For  this  type  of  solution,  the  optimum  focus  would  occur  for 
minimum 


<CF3+CB3 


2 


2 

The  problem  of  minimizing  C„,  with  all  the  other  modes  perfectly  cor- 

“  2 

rected  is  formally  similar  to  that  of  minimizing  oRi  with  only  aF^ 
to  vary.  The  equation  for  zero  derivative  is 


0  = 


(CF3  +  CB3) 


3CF3  ,^B3^CF3  ,3CF33\ 


3N 


3  CF3  3NK 


3Vh  3NK 


The  term  9cB3/3Cp3  can  be  expressed  as 


3  C 


B3 


3  C 


=  C 


F3 


3oE 

3ot 


Where  C  is  the  ratio  of  the  blooming  focus  dispersion  to  the  total 
blooming  dispersion.  Carrying  out  the  differentiation,  we  have 


3  C 


B3 


3  C 


F3 


C°F  °B 
2  2 
°L"°W 


(4/a3)  CF3  CB3 


W 


Which,  being  linear  in  CF3,  leads  to  the  analog  for  the  optimization 
of  opi  : 
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Again,  we  see  the  compensation  cannot  be  perfect,  and  we  see  for  strong 
enough  blooming,  the  optimum  compensation  may  even  be  in  the  same  sense 
as  the  blooming  itself. 

For  future  reference,  note  that  the  optimum  C°3  corresponds  to 
an  optimum  focus  dispersion  0°,  an  optimum  range  ratio  N°,  or  an 
optimum  focal  range  Z^: 

0°  =  Cp3  /2/a 3 

N°  =  l-a°Zt//Om»Rm 

Z°f  =  VNR 

The  optimum  peak  intensity  corresponding  to  C°3  can  now  be  com¬ 
puted.  The  resulting  formula  is  important  because  it  represents  a 
fundamental  limit  and  so  allows  a  prior  judgment  of  the  value  of 
adaptive  optics,  without  regard  to  particular  implementations.  We  begin 
by  evaluating 


Then,  we  form 

4  ■  ■  IC?3*CB3|24 
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and 


and  evaluate 


=  2o.  o. 

Rx  Ry 


PT 


4  >4(1> 


4> 


3 . 2  Limits  of  Thermal-Blooming  Compensation 

Section  3.1  developed  the  fundamental  equations  describing  the 
limits  on  improvement  to  peak  intensity  by  phase  compensation  of  thermal 
blooming,  and  this  section  explores  the  numerical  consequences  of  those 
equations. 


We  begin  by  assuming  the  modes  higher  than  focus  are  fully  cor¬ 
rected.  As  an  example  of  what  this  step  by  itself  may  accomplish,  con¬ 
sider  the  blooming  components  on.  cited  in  Section  2  as  typical  for  a 

til 

uniform  round  beam: 


°B2 

°B3  = 

°B4  = 

If  compensations  are  applied  such  that  there  is  no  net  refocus 


and  such  that  all  modes  higher  than  focus  are  compensated, 

°Ri  =  °Fi  +  0Bi  =  Constant 


°B//I 

CB//2 
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then,  the  constant  must  be 


and  the  compensations  must  be 


The  individual  remaining  are  equivalent  to  an  overall 


=  0.85  o0 


This  is  the  same  thing  as  the  focus  component  of  the  original  set  of 
oD. 's.  The  factor  of  0.85  may  seem  surprisingly  large  for  a  focus 
content  of  thermal  blooming,  due  to  the  fact  that  it  has  been  so 
customary  to  think  in  terms  of  Zernike  coefficients  rather  than  actual 
dispersions.  Indeed,  even  though  the  focus  Zernike  coefficient  may  be 
unimpressive  compared  to  other  modes,  the  relationship 


and  the  numerical  values  of  the  or  on  the  order  of 


°i  "  [°Dk(/Pirdr)1/2] 
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make  focus  the  important  contributor  in  terms  of  dispersion.  Its 
fractional  contribution  to  the  overall  blooming  area  oB  is  (0.85)  =0.73. 

Consequently,  correction  of  modes  higher  than  focus  reduces  the 
deterministic  ray-bending  dispersion  of  the  beam  by  a  factor  of  only  0.73. 

Now,  let  us  proceed  to  apply  refocus  to  achieve  whatever  further 
correction  may  be  possible.  The  preceding  section  provides  a  formula 
for  the  optimal  combination  of  applied  focus  and  blooming  focus  Zernike 
coefficient,  which  when  taken  to  first  order  in  Cg3  leads  to  the  implica¬ 
tion  that  optimal  focus  correction  will  reduce  the  residual  beam  dis¬ 
persion  due  to  deterministic  ray  bending  by  roughly  a  factor  of 


(C_B3  3CB3Y 

yCF3  3CF3 J 


This  can  be  recast  simply  in  terms  of  dispersions  as 


and  is  approximately  equal  to 


/° .73  oj 

I  2  2 

\  L  "  °W 


The  overall  reduction  in  due  to  correction  of  higher  modes, 
and  due  to  optimal  correction  of  focus,  is  thus  approximately 


0.73  (0.73  o2b)2 


(0L-°W) 


In  terms  of  peak  intensity,  the  overall  improvement  factor  is 


2  2 
,  (0.73  Op)3 

°l(1)  +— 2 — hr 
(0L  -  V 
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or 


P  2  2  2  2 

(g£d)  +  q  b)(VV 

a3(l)  (a^-  o^)5  +  (0.73  a*)3 

If  the  blooming  is  small  in  comparison  to  diffraction,  improvement  is 

neither  possible  nor  necessary.  But  if  blooming  dominates,  the  improve- 

2 

ment  ratio  is  very  definitely  limited.  Even  if  we  allow  that  the  oT 

2  L 

includes  not  only  the  random  effects  represented  by  oT(l),  but  also 

the  deterministic  effects  of  the  applied  oF^»  and  even  if  the  latter 

are  assumed  to  roughly  equal  minus  the  o_.,  the  improvement  ratio  cannot 
3  Bl 

exceed  l/(0.73)  ~  2.6.  This  result  confirms  the  limitation  that  other 

researchers  have  previously  encountered  in  simulations  and  tried,  more 
or  less  fruitlessly,  to  improve  upon. 

One  technique  to  improve  the  result  a  little  bit  is  to  change  the 
problem  geometry  in  a  way  that  will  change  the  spreading  in  the  y  di¬ 
mension  from  the  nearly  uniform  situation  represented  by  Og^  =  cr^  = 

Og//2"  to  one  with  more  disparity — say,  by  introducing  a  nonuniform  or 

noncircular  beam.  If  such  manipulation  could  get  to  the  limit 

oB3  ~  0,  o04  ~  Og,  then  correction  of  modes  higher  than  focus  would  result 

in 


Accordingly,  the  factor  of  0.73  changes  to  0.5,  and  the  upper-bound 
correction  ratio  changes  to  (1/0. 5) 3  =  8.  The  factor  of  8  seems  much 
more  impressive,  except  when  viewed  in  the  perspective  that  it  would 
occur  only  in  situations  that  were  catastrophically  bloomed  in  the 
first  place. 

We  can  say  in  general,  that  if  C  is  the  ratio  of  blooming  focus 
dispersion  to  total  blooming  dispersion,  then  the  upper  bound  of 

3 

blooming  compensation  is  1/C  .  This  result  constitutes  a  very  strong 
statement  because  it  is  not  at  all  scenario  or  hardware  implementation 
dependent.  The  only  ingredients  of  the  analysis  are  the  basic  form  of 
the  algebraic  model  and  the  basic  geometry  of  the  ray-propagation  process. 
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3.3  Implications  to  Linear  Compensation 


It  seems  intuitively  obvious  that  maximum  peak  intensity  ought 
to  result  if  beam  spread  due  to  turbulence  and  jitter  is  made  as  small 
as  possible.  This  supposition  turns  out  to  be  naive,  however,  given 
the  limitations  on  thermal-blooming  compensation  cited  in  the  last 
section.  The  problem  is  that  the  total  beam  degradation  includes  both 
the  residual  "high-frequency"  beam  spread 


2  2  .  2 
aH  =  °T  +  0J 


and  the  optimal  ray-bending  dispersion 


(C  o|>3 


which  depends  upon  through  at  -  o £  and  o^.  The  is  a  linear 

j  n  Jj  W  d  H 

contributor  to  -  o-. 


2  2 
°L'aW 


2  2  2 
°D  +  °H  +  °B 


2  2 

Here  a  is  the  diffraction  term,  and  oa  approximates  the  aberration 

u  B  2  0 

compensations  forced  at  the  aperture.  Then,  the  a„  also  enters  at 

H  B 

through  the  definition 


»B  '  Ci<°L-°X 


Not  only  is  it  there  explicitly  in  -  a*,  but  also  in  the  Yh  through 
the  parameters  and  noted  in  Section  2.2. 

Optimal  linear  compensation  has  to  mean  minimization  of  the 
combination : 


2  ,  2 
°H  +  °R0 
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This  may  occur  for  nonzero  ou 

n 


if  the  derivative 


3  o 


1  + 


RO 


9o; 


is  zero  for  a  nonzero 
2 


2  .  2 
o„  (or  negative  for  zero  o„) 

ri  n 


oRQ  can  be  simplified  to 


The  expression  for 


°R0  '  ‘CCi'h>3  <0L-«w' 


so  the  required  derivative  is 


l  +  o 


2 

RO 


where 


3Vnq  +  3Vnq 

>0H 

3Nq3o2  3N’3a2 

will  be  a  negative  number,  possibly  allowing  the  required  derivative 

2 

to  be  zero  for  nonzero  o„. 

H 

To  examine  the  situation  more  quantitatively,  we  required  the 
derivatives 


2  2  2 
2<m2  o2  +  o2) 


3NQ  =  No 

'  2  (oj-oj) 


and  the  corresponding  derivatives  of  These  are  difficult  to  evaluate 

unless  we  introduce  some  approximations.  If  the  Fresnel  number  Np  is 
large  and  the  range  ratio  NR  is  close  to  1,  we  have 


exp 


and 


NN_ 
D  F 


nq<ns  +  1) 


The  required  derivative  can  be  evaluated  in  terms  of  the  above, 
and  for  large  N  ,  reduces  approximately  to 

r 

1  +  (C  tj) 3  +  1] 


In  any  situation  where  the  above  is  negative,  perfect  linear  correction 
is  not  optimal.  The  existence  of  such  situations  is  a  surprising 
result.  But  it  seems  to  suggest  an  explanation  for  what  has  previously 
been  observed  in  simulations}  that  the  presence  of  even  a  small 
amount  of  thermal  blooming  can  make  turbulence  correction  unstable. 

The  above  formula  depends  only  on  heating  phase  and  not  on  actual 
blooming  dispersion  a3.  Even  if  the  latter  seems  negligibly  small,  the 
former  may  be  large  enough  to  make  full  linear  correction  suboptimal. 

Naturally,  one  wants  a  generally  applicable  rule  for  defining 
optimal  turbulence  and  jitter  correction,  assuming  optimal  thermal- 
blooming  compensation  is  available.  The  rule  follows  from  the  fact  that 
we  would  like  to  end  up  with  the  heating  phase  no  larger  than 


hO 


(^-D1/3 


CC^ 


1/a 
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2 

after  driving  the  high-frequency  beam  spread  to  an  optimal  value  crH0» 
and  applying  the  optimal  blooming  correction.  Let  us  suppose  that 
the  sensitivity  of  heating  phase  is  mainly  in  the  factor,  and  that 
N'  without  corrections  is  proportional  to 


2  2 
aD  +  0H 


and  at  optimal  correction  to 


°D  +  °H  0  +  °B 


2 

where  aD  approximates  the  aperture  corrections  for  blooming.  This  is 
a  2 

not  explicitly  known  until  is  already  specified,  so  we  proceed  by 

2H0  (4) 

further  approximating  the  ofi  in  terms  of  the  algebraic  model.  This 
makes  N'  corrected  proportional  to 


m 


2 


This  means  high-frequency  beam  spread  should  be  made  no  smaller 
than 


2  ,  2  2X  h  2  2 

°H0  =  °D  +  °H  ^  -  m  aD 

if  that  is  a  positive  number. 

3.4  Implications  for  Power  Throttling 

The  concept  of  varying  laser  power  for  optimization  has  been 

discussed  elsewhere ^  and  shown  to  be  exceedingly  sensible.  It  seems 

to  have  a  clear  application  in  regard  to  adaptive  optics.  In  Sections 

2 

3.2  and  3.3,  we  have  optimized  the  ray-bending  dispersion  o  for  fixed 

2  K 

high-frequency  beam  spread  oH  with  P  held  fixed,  and  now  finally  we 
can  consider  optimizing  P. 
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The  ultimate  objective  is  to  optimize  the  peak  intensity 


I 

P 


2  2  2  2 
lrZt(°D  +  0W  +  0HO  + 


'RO 


) 


2  2 

given  the  functional  dependence  of  o  n  and  oun  upon  P.  Using  the 

(4)  KU  MU  2 

algebraic  model  parameter  m  to  estimate  the  oD  that  would  not  be 

D 

known  prior  to  specifying  P,  we  have 


(CCB’hO|3 


(m 


°D+CJH0) 


where 


( 4+4 ) 


2 


m 


a 


2 

D 


if  that  is  positive,  zero  otherwise,  and  where 


hO 


GM 


1/3' 


c  c  * 


1/a 


Since  is  linear  in  P,  1^  essentially  has  a  denominator  which  is  a 
polynomial  in  P,  possibly  going  up  to  fourth  order.  Analytic  optimiza¬ 
tion  of  the  full  expression  looks  intractable,  but  a  quick  one-dimensional 
computer  search  should  suffice  to  find  a  suitable  operating  point. 
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SECTION  4 


CONVENTIONAL  ADAPTIVE  OPTICS  SCHEMES 


The  fact  is  that  no  existing  adaptive  optics  scheme  fully  realizes 
even  the  modest  corrections  available  for  thermal  blooming,  and  each 
may  actually  go  seriously  wrong.  This  section  shows  what  goes  wrong 
and  how  it  may  be  remedied. 


4 . 1  What  Goes  Wrong 

The  classic  techniques  for  achieving  adaptive  optics  include  open- 
loop,  predictive  compensation,  and  at  least  three  versions  of  closed- 
loop  compensation.  The  first  of  these  is  referred  to  as  "phase  conjugate". 
A  point  source  at  the  receiver  gives  rise  to  a  phase  profile  at  the 
aperture,  which  is  exactly  compensated  by  the  adaptive  mirror  so  that, 
by  reciprocity,  the  output  beam  will  propagate  to  that  point.  The  next 
technique  is  referred  to  as  "multidither".  The  adaptive  mirror  is 
dithered,  either  by  zones  or  by  Zernike  modes,  and  the  radiation  reflected 
from  the  receiver  is  monitored  to  determine  adjustments  in  the  mirror. 
Another  technique  is  referred  to  as  "image  sharpening".  It  is  similar 
to  multidither,  except  that  the  sharpness  of  the  received  image,  rather 
than  the  intensity  of  reflected  radiation,  is  monitored.  The  problems 
with  each  of  the  techniques  are  discussed  in  the  following  paragraphs. 


4.1.1  Predictive  Compensation 

Given  all  the  relevant  variables,  it  is  possible  to  compute  the 

entire  heating-phase  profile  to  be  expected  from  an  uncorrected  beam. 

Normally,  what  is  done  with  this  is  to  set  the  applied  phase  ¥  equal  to 

app 

that  for  a  spherical  wave  converging  on  target,  p^us 

conjugate  of  the  computed  heating-phase  profile,  ^pp5 


T 

app 


-fsph<Bt) 


+  V 


hpp 
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The  static  thermal-blooming  model  indicates  that  one  should  not  do  this. 
Instead,  one  should  go  through  the  following  steps: 

(1)  Compute  optimum  focus,  and  assume  the  optimum  focus  is 
applied. 

(2)  Compute  the  phase  that  would  then  result,  and  assume  modes 
higher  than  focus  are  fully  corrected. 

The  difference  between  these  two  strategies  is  that  the  conven¬ 
tional  one  always  applies  more  focus  than  ¥Sph(Zt)»  whereas  the  newly 
suggested  one  can  apply  less  focus  than  ^gp^^^.)  in  case  of  heavy 
blooming. 

It  is  noteworthy  that  one  of  the  major  motivations  for  closed- 
loop  adaptive  optics  has  been  the  lack  of  any  proof  that  open-loop 
correction  is  optimal.  The  previous  argument  indicates  that  the 
conventional  open-loop  correction  is  indeed  not  optimal,  but  the 
analysis  also  suggests  how  to  optimize  it  without  having  to  resort  to 
a  closed-loop  fix — such  as  dithering  the  amplitude  of  ¥ 


4.1.2  Phase  Conjugate  Adaptive  Optics 


Phase  conjugate  adaptive  optics  is  normally  implemented  simply 


by  driving  ¥  to  the  conjugate  sensed  phase  T 

app  3  sen 


But  the  static 


thermal -blooming  analysis  provides  some  additional  information  that 
can  be  used.  We  know  the  optimal  focus  correction  is  and 

is  clear  that  if  the  optimal  correction  were  applied,  the  conjugate 
sensed  phase  would  be  fgen  =  ^gp^CZ^)  +  ^hpp •  Let  ^hpf  be  tbe  focus 
part  of  ^hpp"  We  can  define  the  ratio  of  focus  phases 


’sph(z?» 


Vh'V  ♦Vf11# 


Apparently,  optimal  correction  cannot  be  found  by  driving  t  to  f 

app  sen 

Instead  it  requires  driving  only  the  higher  modes  to  such  an  equality. 

For  focus  parts,  one  should  drive  ¥  to  N°.  4* 

app  R  sen 


The  previous  observations  can  help  to  explain  and  remedy  some 
mysterious  results  that  have  been  thought  to  occur  with  conventional 
phase  conjugate  adaptive  optics.  Figure  4  shows  three  major  categories 
of  what  a  phase  conjugate  system  can  do,  according  to  the  nonlinear 


3 


Figure  4.  Conventional  phase  conjugate  results. 

lens  models  of  Hermann ^  and  the  computer  simulation  of  Brown,  et  al!2^ 
These  three  possible  phase  conjugate  results  can  be  identified  with 
three  possible  circumstances  that  can  occur  according  to  the  static 
thermal-blooming  analysis: 

(1)  N°  >  1,  N°  >  1;  the  optimum  focus  H*  >  ¥  is  almost 

i\  x\  cipp  sen 

reached,  but  not  quite.  This  situation  can  be  called 
light  blooming. 

(2)  N°  >  1,  N°,  <1;  the  system  walks  through  the  optimum 

focus  correction,  which  occurs  for  Y  <  t  This 

app  sen 

situation  can  be  called  moderate  blooming. 

(3)  N°  <1,  N^,  <  1;  the  system  goes  wrong  way,  focusing  in 

when  the  optimum  is  behind  the  target.  This  situation 
constitutes  heavy  blooming. 
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The  previous  analysis  means  that  the  optima  in  the  three  phase 
conjugate  situations  are  as  indicated  in  Figure  5.  In  no  case  is  the 
system  achieving  and  holding  an  optimum! 


Figure  5. 


Conventional  results  in  relation  to  optimum. 


It  is  noteworthy  that  the  terms  light,  moderate,  and  heavy  are 
defined  by  the  values  of  N°  and  N°,,  and  not  by  the  Strehl  ratio 

K  K  O  O 

I  /I  ,  where  u  denotes  unbloomed.  Because  N°  and  N„.  depend  on  a 
p  pu  R  R 

single  phase,  whereas  I  /I  depends  on  a  whole  phase  profile,  it  is 

1.  ir 

possible  for  heavy  blooming  to  have  a  large  Strehl  ratio  I  /I 

p  pu 

The  behavior  already  described  has  not  previously  been  subject  to 
remedy.  Instead,  the  poor  behavior  has  been  attributed  to  uncorrect- 
ability  of  blooming  far  from  the  aperture.  Indeed  such  blooming  is 
somewhat  uncorrectable,  but  this  model  indicates  that  the  uncorrect- 
ability  only  determines  what  the  optimum  peak  intensity  is,  and  not 
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why  it  is  not  achieved  and  held.  Also,  poor  behavior  has  previously 

been  attributed  to  the  use  of  a  linear  control  law  in  a  nonlinear 

system.  Indeed,  the  linear  control  law  f  =  f  _  does  not  work, 

app  sen  - 

but  only  because  it  has  been  applied  to  all  modes  equally. 

4.1.3  Multidither  Adaptive  Optics 

The  intent  of  the  multidither  scheme  is  to  drive  N  to  N°  by 

R  K 

observing  the  time  average  value  of  the  product  of  commanded  phase 
deviation  and  observed  intensity  deviation.  In  principle,  the  system 
should  work;  but  there  are  several  practical  problems.  The  system 
tends  to  have  difficulty  if  started  very  far  from  the  solution  it 
seeks,  and  tends  to  provide  inadequate  information  concerning  the 
size  of  steps  required.  There  are  other  problems  as  well,  but  those 
mentioned  are  potentially  reducible  on  the  basis  of  knowledge  gained 
from  the  static  thermal-blooming  analysis.  Figure  6  illustrates  the 
problems . 

VANISHING  STEPS 


Figure  6.  Illustration  of  some  multidither  problems. 

The  analysis  is  helpful  because  knowledge  of  N°  allows  the  system 
to  be  initialized  predictively  near  the  optimum,  and  having  a  model 
for  Ip(N  )  means  step  size  can  be  properly  chosen.  To  control  the  step 

R 
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size,  it  is  useful  to  measure  the  curvature  as  well  as  the  slope  of  I 

P 

versus  Vcam  by  using  the  harmonics  of  the  return  signal.  Then,  the  slope 
and  curvature  can  be  used  to  define  a  parabola.  If  the  parabola  has  a 
peak,  the  step  should  go  directly  to  it.  If  not,  the  ratio  of  the  first 
and  second  derivatives  may  provide  some  guidance. 


4.1.4  Image-Sharpening  Adaptive  Optics 


Finally,  let  us  consider  adaptive  optics  by  the  image-sharpening 
technique.  This  technique  is  subject  to  the  same  suboptimal  equilibrium 
as  the  phase  conjugate  technique,  plus  the  same  stalls  and  overshoots 
as  the  multidither  technique.  All  of  the  same  remedies  are  required. 


The  hope  for  remedying  any  of  the  conventional  adaptive  optics 
approaches  in  regard  to  thermal  blooming  resides  in  techniques  for 
determining,  from  available  data,  what  the  optimal  focus  correction  really 
is.  If  it  can  be  determined,  then  predictive  compensation  can  be  done 
properly,  phase  conjugate  compensation  can  be  modified  to  stabilize 
at  the  optimum,  multidither  can  be  properly  initialized,  and  image 
sharpening  can  be  both  initialized  and  kept  from  going  to  a  suboptimal 
phase  conjugate  solution. 


The  key  is  to  try  to  identify  sufficient  data  to  fill  in  all  the 
ingredients  of  the  equation  for  optimal  focus  correction 


Let  us  suppose  a  thermal-blooming  phase  profile  is  available, 
either  predicted  or  extracted  from  a  return  wave  by  an  interferometer 
or  an  image-sharpening  mirror.  That  is,  let  us  assume  we  have  access 
to  characteristics  like  blooming  Zernike  coefficients  and  phase  perturba¬ 
tion  at  the  center  of  the  aperture. 

Let  us  suppose  also  that  some  basic  system  parameters  are  available, 
including : 
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_ . iC— .  i - jJ?.'- 


2  ,  2 
°T  +  0J 


aperture  radius 
wavelength 
target  range 
beam  quality 

high-frequency  tilt  statistics 
low-frequency  tilt  statistics 


(4) 

The  algebraic  propagation  model  provides  some  key  numbers 


C ' 

a 

m 

m ' 

m" 

Infinite  Gaussian 

0.010590 

1.6419 

2.0 

0.6366 

1 

Truncated  Gaussian 

0.028727 

1.3715 

4.0 

0.9166 

0.8893 

Uniform  Beam 

0.014264 

1.1777 

1.0 

0.9202 

1.124 

and  some  key  equations 


=  2tt  R„/X  Z, 


=  0.5(X/itRm)1 


=  0. 5 (m'MX/2R  )' 

m 


k  =  2tt/X 
2  2  . 


2  2  2  2  2  2 
=  oD +  aT + Oj +  aw +  oB  (oB  =  corrections) 


=  (m 


2  cl.  cl,  A1'2, 


Section  2  of  this  report  provides  the  definition  of  conversion 
factors 


«i  ■  %  k(/p?tdr)l/2 


37 


which  allow  evaluation  of  the  blooming  dispersions 


°B  "  °D  £  1CBi^i>2 

and,  consequently,  the  fraction  of  total  blooming  area  due  to  focus 
aberration 


C 


CTD(CB3/a3r 


(which  we 
The 

which  can 


already  know  from  Section  3  to  be  on  the  order  of  0.73). 

(4) 

algebraic  propagation  model  depends  upon  a  heating  phase 
be  defined  from  the  data  as 


’h  ‘  [°B/CB  l^-4>]1/a 

From  the  definition  of  4^  shown  in  Section  2,  it  must  also  be  rather 
similar  to  (N^)  1  times  the  phase  perturbation  existing  at  the  center 
of  the  aperture,  so  that  pieces  of  data  may  provide  an  alternative 
basis  for  evaluation. 

To  fill  in  the  equation  for  the  optimal  focus  correction,  we  also 
need  a  derivative  of  this  phase:  3  4^/3  N  .  Using  any  of  the  analytic 
approximations  for  4^,  3  4*^/3  N  can  be  approximated  analytically.  For 
example,  suppose  we  consider  4'h  in  the  limit  of  large  Fresnel  number 


’h  *  lND»FVNQ|,,StV|e,t',(-T1) 


In 


2(NS  +  1) 


(1  -  V  * 


[(1  -  »R 


>*♦  <2N„  W2] 
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Then,  3  4^/31^  is  obtained  as  follows.  First,  we  note  that  the  factors 
ND,  Np,  and  NR  are  each  individual  functions  of  NR  for  a  given  Zt,  but 
the  dependencies  cancel  leaving  the  product  NQ  NR  NR  independent  of  NR. 
Thus,  3  4*^/3  Nr  has  only  three  contributors 


ns  +  nr  + 


4* 

Th  3  In 
In  3Nr 
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where 
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One  may  thus  conclude  that: 

(1)  Predictive  compensation  is  readily  achievable  if  the  blooming 
phase  profile  can  be  adequately  predicted. 

(2)  Phase  conjugate  compensation  can  be  suitably  modified  if 
blooming  can  be  separated  from  turbulence. 

(3)  Multidither  compensation  can  be  helpfully  initialized  if 
the  blooming  phase  profile  can  be  predicted. 

(4)  Image-sharpening  compensation  can  be  helpfully  initialized 
and  suitably  modified  if  blooming  can  be  predicted  and  then 
separated  from  turbulence. 

So  far  as  closed-loop  adaptive  optics  is  concerned,  the  key 
problem  remaining  is  separation  of  blooming  from  turbulence,  which  has 
to  be  simultaneously  sensed.  Dynamic  behavior  is  the  only  basis  for 
such  a  separation. 
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SECTION  5 


DYNAMICS 


The  existence  of  temporal  dynamics  is  both  an  assistance  in  the 
design  of  an  operational  adaptive  optics  system  and  a  confusing  factor 
in  the  prediction  of  its  performance.  This  section  discusses  both 
aspects  of  the  situation. 

5 . 1  Discrimination  Between  Phenomena 

Section  4  presented  a  complete  algorithm  for  optimal  thermal 
blooming  correction,  leaving  as  a  key  requirement  the  discrimination 
between  turbulence  and  thermal  blooming.  Dynamic  behavior  provides 
the  basis  for  such  discrimination. 

Taking  a  first  cut  at  the  discrimination  problem,  one  may  note 
that  there  tends  to  be  a  difference  in  frequency  content  between  the 
two  phenomena.  Each  can  be  thought  of  as  a  set  of  phase  screens  moving 
past  the  aperture  due  to  wind  and  slew  convection,  while  also  evolving 
slowly  due  to  conduction  as  they  go.  The  blooming  phase  screens  are 
zero  everywhere  except  where  deposited  heat  has  caused  density  perturba¬ 
tion;  whereas,  the  turbulence  phase  screens  continue  indefinitely  on  a 
spatial  scale  of  p,  the  coherence  length.  This  means  that,  in  a  condi¬ 
tion  of  high  slew  and/or  small  coherence  length,  changes  in  turbulence 
profile  will  be  on  a  shorter  time  scale  than  those  in  the  booming  profile. 

Now,  suppose  we  have  to  face  the  more  difficult  situation  where 
the  time  scales  for  blooming  and  turbulence  indicate  temporal  over¬ 
lapping  of  the  two  phenomena.  There  is  still  a  basis  for  discrimination 
in  the  spatial  peculiarity  of  the  blooming  profile.  The  main  peculiarity 
is  that,  either  by  convection  or  buoyancy,  the  blooming  tends  to  move 
across  the  aperture,  being  zero  on  one  side  and  nonzero  on  the  other. 
Normally,  the  zero  edge  is  leading  into  the  wind,  and  in  any  event  it 
can  be  located  on  the  basis  of  the  time  average  squared  phase  amplitude 
at  the  edge. 
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The  phase  profile  history  along  the  edge  where  blooming  should  be 
zero  can  provide  a  good  estimate  of  turbulence  everywhere  later  on. 

The  trick  is  to  break  the  observed  edge  profile  into  separate  phase 
sheets  corresponding  to  different  transverse  velocities  and  reassemble 
them  with  appropriate  time  shifts  for  other  points  in  the  aperture.  For 
convenience,  suppose  this  edge  is  leading  into  the  wind  and  the  wind  is 
along  x.  Let  -x  be  the  x  coordinate  of  a  point  on  the  edge.  The  phase 
here  is  H*  (~xm,  y,  t).  The  differential  phase  screen  moving  at  velocity 
v  can  be  identified  as  the  frequency  component  of  t  at  frequency  f  =  v/p. 
The  turbulence  phase  at  x  can  be  estimated  by  multiplying  the  frequency 
component  of  ¥(-xm,  y,  t)  at  frequency  f  by  the  phase  shift  factor 
ei2Trf\  where  t  =  (-xm+x)/v,  and  resynthesizing. 

5.2  Modeling  Expected  Behavior 

All  of  the  conventional  adaptive  optics  algorithms  and  the  new 
ones  here,  suggested  for  optimal  thermal-blooming  correction  and  for 
turbulence-blooming  discrimination,  are  apt  to  exhibit  complex  behavior 
in  a  dynamic  situation.  In  order  to  anticipate  what  may  happen,  it  is 
necessary  to  have  some  models. 

Let  us  begin  with  the  dynamics  of  thermal  blooming.  Appropriate 
models  can  be  based  upon  the  algebraic  model  developed  so  far,  but  the 
temporal  variation  requires  some  modification  to  the  quantities  in¬ 
volved.  An  extended  model  incorporating  dynamics  must  take  account  of 
the  fact  that  the  atmosphere  has  some  memory  due  to  the  finite  time 
required  for  heated  air  to  pass  out  of  the  propagation  path,  and  it 
must  take  account  of  the  fact  that  the  corrections  actually  applied 
may  lag  commands  somewhat  because  they  are  produced  by  some  type  of 
servo  mechanism. 

Atmospheric  memory  will  appear  in  an  effective  heating  phase 
representing  some  type  of  time  average  over  the  recent  past.  This 
effective  phase  must  impact  the  instantaneous  peak  intensity  Ip(t) 
directly,  and  also  indirectly  through  the  instantaneous  thermal-blooming 
beam  spread  parameter  m.  Figure  7  is  a  functional  block  diagram  of  the 
relationships  in  the  dynamic  system. 

The  actual  computations  involved  in  most  of  the  blocks  have 
already  been  described  in  Section  2.  Only  ?h  and  m(t)  remain  to  be 
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LASER. 

SCENARIO. 

ATMOSPHERE 

VARIABLES 


Figure  7.  Dynamic  model  of  propagation  process. 

defined.  Let  the  be  the  time  average  over  the  recent  past,  with 
the  beam  profile  at  the  aperture  as  a  weighting  function: 


R 


Let  the  m(t)  be  contrived  to  have  reasonable  limiting  values: 

m(t)  =  1  before  blooming  develops 

m(t)  =  m  at  steady  state 

m(t)  >  m  if  blooming  has  been  heavier 

in  the  past. 

A  simple  form  for  m(t)  is 

m (t)  =  1  +  (m  -  1) aB/oB 

Consider  now  the  imposition  of  any  of  the  standard  adaptive  optics 
algorithms  on  the  system  described  by  Figure  7.  In  the  case  of 
predictive  compensation,  the  aFi  simply  correspond  to  constants  equal 
to  minus  the  values  the  would  have  in  uncorrected  steady  state. 
Imposition  of  any  of  the  closed-loop  schemes,  however,  is  more  compli¬ 
cated,  requiring  an  additional  loop  in  the  system. 
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Figure  8  illustrates  a  phase  conjugate  system.  For  this  system, 
the  adaptive  corrections  correspond  to  whatever  phase  is  sensed,  which 
in  turn  depends  on  the  atmospheric  heating.  To  minimize  dealing  with 
whole  profiles  instead  of  simple  numbers,  let  us  adopt  the  following 
simple  model.  For  the  sensed  phase,  suppose 


sen 


V^t^h 


This  would  correspond  to  a  focus  correction 


R  4\, 

/7T— C  m  _ *L 

op  -  -/0.5  z  y 


Jt  sph (Zfc) 


Partitioning  into  individual  Op^  in  the  usual  way. 


°F1  =  ° 


°F2  ~  °F 


°F3  *  °F4  ^  °F//2 


LASER, 

SCENARIO, 

ATMOSPHERE 

VARIABLES 


Figure  8.  Dynamics  of  phase  conjugate  compensation. 


Figure  9  shows  a  multidither  system.  For  this  system,  oscillating 
perturbation  corrections  Aa_.  are  observed  via  a  glint  return,  and  some 

r  1 

update  to  a„.  is  rendered  on  the  basis  of  how  they  compare.  The  equations 
FI 

can  be  modeled  as 


c  <  Aof.  AIp  > 


I  f  I  Ee(t_T)  E-dT 

P  J  P 

T  =  0 
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Figure  9.  Dynamics  of  multidither  system. 

Figure  10  shows  an  image-sharpening  system.  Here,  oscillating 
perturbation  corrections  AoFi  are  entered,  and  the  resulting  target 
image  is  monitored  for  oscillating  sharpness  S,  with  some  update  to 
oFi  being  based  upon  the  comparison  of  the  two.  (Monitoring  sharpness 
is  essentially  equivalent  to  monitoring  phase.)  The  equations  are 
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Finally,  consider  the  impact  of  servo  lag  on  all  the  systems. 

The  a  actually  applied  are  a  result  from  a  servo  mechanism  trying  to 
implement  the  commands.  Therefore,  residual  discrepancies  will  always 
exist,  usually  in  the  form  of  lags.  An  example  is  illustrated  in  Fig¬ 
ure  11:  a  type  zero  servo  drives  a  to  o_ .  .  For  this  example, 

Fi  app  Fi  com  r 

the  frequency  domain  equations  are 


°Fi  app 


|  °Fi  com^  °Fi  app 


..  X 
s  +  y 


Figure  11. 


Example: 
driving  a 


type-zero  servo 

to  o_ .  . 

Fi  com  Fi  app 


and  the  time-domain  equations  are 


°Fi  app(t) 


°Fiapp<0) 


h  J  °Fi  com 1 
0 


r  j.  app 


The  equations  reveal  the  lag:  whatever  aFi  is  commanded  after  t=0  is 

not  fully  applied  for  a  while.  More  generally,  the  main  point  of  any 

such  model  is  that  o„.  (t)  laqs  a„.  (t)  because  an  integration  is 

ri  app  Fi  com 

involved. 


5.3  Implications 

Inclusion  of  temporal  variation  adds  some  nuances  to  the  analyses 
of  what  each  of  the  adaptive  optics  systems  will  do,  and  how  to  improve 
the  behavior.  Generally,  the  dynamics  will  cause  overshoots  and  oscilla¬ 
tions,  and  sometimes  even  wrong  commands. 
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Consider  the  case  of  predictive  compensation.  When  the  system 
is  first  turned  on,  peak  intensity  is  lower  than  it  will  become,  so 
blooming  is  less  than  it  will  become.  Therefore,  the  peak  intensity 
will  jump  initially  (on  the  time  scale  of  the  servo)  to  some  value,  and 
and  then  it  will  either  climb  or  decline  slowly  (with  the  atmospheric  time 
scale)  toward  equilibrium.  In  either  case,  it  might  overshoot  and 
oscillate,  as  shown  in  Figure  12. 


t 


Figure  12.  Oscillations  superposed  on 
basic  predictive  behavior. 

Next,  consider  the  case  of  phase  conjugate  adaptive  optics.  For 

moderate  and  heavy  blooming,  lp(t)  is  declining  as  it  approaches  its 

equilibrium  value.  Because  Ysen  lags  Ip,  '*'sen  is  then  larger  than 

at  equilibrium;  and  because  y  lags  it  is  then  even  larger 

app  3  sen' 

than  the  equilibrium  1*'sen-  Consequently,  the  system  overshoots,  and 
oscillations  follow,  as  illustrated  in  Figure  13. 


Figure  13.  Oscillations  superposed  on  basic 
phase  conjugate  behavior. 

Now,  consider  multidither  adaptive  optics.  The  scheme  uses  the 

expected  value  <AI  Af  >  to  drive  N„  to  N°.  But  servo  lag  and 
p  com  R  R 

atmospheric  integration  time  can  significantly  change  the  phase  between 
Alp  and  AfCQm,  possibly  giving  the  wrong  sign!  This  will  at  least 
delay,  and  may  even  prevent,  convergence. 

To  encourage  proper  functioning  of  a  multidither  scheme,  it  is 

helpful  to  try  to  overcome  the  lags.  One  approach  is  to  dither  as 

slowly  as  practical.  This  approach  is  of  limited  value  because  it 

compromises  the  bandwidth  capability  of  the  system.  Another  approach 

is  to  attempt  to  measure  the  phase  of  AI  with  respect  to  AH'  .  This 

p  com 

can  be  difficult,  and  2n  ambiguities  are  a  problem.  Further  assistance 
can  be  gained  by  using  at  least  two  dither  frequencies  per  mode  and 
extrapolating  to  zero  dither  frequency. 

Next,  consider  image-sharpening  adaptive  optics.  It  has  the  same 
problems  due  to  phasing  between  perturbation  introduced  and  perturbation 
observed  as  does  multidither,  and  perhaps  worse,  because  it  dithers  a 
low-power  mirror  which  has  better  bandwidth  capability  than  the  high- 
power  mirror  used  by  multidither.  But  the  problems  shared  with  multi¬ 
dither  can  be  overcome  similarly. 
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Now,  consider  systems  augmented  with  the  algorithm  for  optimal 
thermal-blooming  correction.  Oscillations  will  still  occur  as  in  the 
conventional  systems.  However,  the  oscillations  are  informative:  their 
frequency  doubles  if  the  algorithm  is  optimal  instead  of  conventional, 
because  then  the  equilibrium  about  which  they  oscillate  is  a  maximum. 

Finally,  consider  a  situation  where  turbulence  and  blooming  are 
present  simultaneously.  Failure  to  discriminate  could  be  serious 
because 


| under 
I  over 
results  in 
|  over 
I  under 


I 

I 

I 

1 


estimation  of  blooming 


correction  of  both  phenomena 


and  hence  suboptimal  results. 
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SECTION  6 


CONCLUSIONS  AND  RECOMMENDATIONS 


This  report  argues  that  the  nonlinear  phenomenon  of  thermal 
blooming  presents  a  serious  problem  for  continuous-wave  adaptive  optics. 
The  potential  for  blooming  correction  is  fundamentally  limited.  If 
adaptive  optics  is  to  be  used  at  all  for  this  type  of  system,  the 
motivation  must  lie  in  the  potential  for  correction  of  simpler  linear 
effects.  But  if  adaptive  optics  is  used  for  those  purposes,  blooming 
has  to  be  estimated  and  treated  intelligently,  because  naive  approaches 
will  result  in  catastrophic  failure  to  do  even  the  simple  task  of 
linear  correction,  just  because  the  blooming  is  present.  If  this  is 
true,  then  it  is  obvious  that  adaptive  optics  has  not  yet  been  imple¬ 
mented  with  a  proper  understanding.  Forging  ahead  without  that  under¬ 
standing  would  be  wasteful. 

Verifying  the  analyses  presented  here  requires  that  a  number  of 
key  assertions  and  conclusions  be  confirmed  by  further  analysis,  by 
simulation,  or  by  experiment.  It  is  important  to  break  the  verification 
problem  into  as  many  independent  questions  as  possible,  in  order  that 
each  one  be  as  simple  and  straightforward  as  possible.  The  following 
partitioning  of  the  problem  is  recommended: 

(1)  In  an  unbloomed  beam  with  aberrations  of  focus,  astigmatism, 
primary  coma,  and  spherical  aberration,  the  Zernike 
coefficients  imply  beam-spread  parameters  that  accurately 
predict  peak-intensity  Strehl  ratio  below  beam  breakup. 

(2)  In  an  uncorrected  bloomed  beam  below  pcrit»  the  ratio  of 
defocused  beam  spread  to  total  beam  spread  is  a  function  of 
beam  shape  only. 

(3)  For  a  uniform  round  beam,  the  ratio  of  defocused  beam  spread 
to  total  beam  spread  is  approximately 


°B 
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(4)  Application  of  any  reasonable  blooming  correction  strategy 
has  the  effect  of  increasing  the  value  of  pcr^t»  as  des¬ 
cribed  in  Reference  4. 

(5)  If  arbitrary  Zernike  coefficients  C„.  are  applied  at  the 

r  X 

aperture  and  resulting  blooming  coefficients  CD.  are  sensed, 

til 

the  sum  of  Zernike  coefficients 


CFi  +  CBi 


determine  a  value  of  beam  spread  on  target  that  accurately 
predicts  the  peak-intensity  Strehl  ratio  so  long  as 

P  <Pcrif 

(6)  The  combination  of  uncorrected  focus  and  coma  aberrations 
can  result  in  tilt  aliasing  that  would  in  principle  affect 
hot-spot  tracking. 

(7)  The  optimum  thermal-blooming  correction  strategy  is  to 
fully  correct  modes  other  than  focus,  and  to  optimize 
the  focus. 

(8)  The  optimum  focus  correction  is  not  generally  equal 

to-CB3-  - 

(9)  It  is  possible  for  the  optimum  to  have  the  same  sign 

as  CB3. 

(10)  The  numerical  value  of  the  optimum  CF3  can  be  estimated 
from  the  algorithm  cited  in  this  report. 

(11)  The  limits  on  blooming  correction  are  a  function  of  beam 
shape  only,  and  amount  to  a  factor  of 


3 


(12)  If  there  are  correctable  linear  dispersions  (jitter  or 
turbulence)  as  well  as  blooming,  in  some  cases,  maximum 
peak  intensity  requires  leaving  some  of  the  linear  spread 
uncorrected. 
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(13)  The  optimum  linear  spread  to  leave  uncorrected  can  be 
estimated  with  the  algorithm  cited  in  this  report. 

(14)  If  blooming  and  linear  effects  are  all  to  be  optimally 
corrected,  there  is  then  an  optimum  power  to  operate 
at,  as  described  in  this  report. 

(15)  Discrimination  between  blooming  and  turbulence  can  be 
accomplished  on  the  basis  of  temporal  frequency  content 
when  slew  is  high  or  coherence  length  is  small. 

(16)  Even  if  blooming  and  turbulence  frequency  spectra  overlap, 
discrimination  can  be  accomplished  in  the  manner  described 
in  this  report. 

The  analyses,  simulations,  and  experiments  aimed  at  verifying 
these  items  need  to  be  designed  in  the  least  burdensome  way  possible. 
It  should  be  possible  to  concentrate  on  geometric  analyses,  limited 
wave-optics  simulations,  and  small-scale  thermal-blooming  cell- 
laboratory-type  experiments. 
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